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ABSTRACT 

We present our second paper describing multi-wave-band, time-resolved spec- 
troscopy of WZ Sge. We analyze the evolution of both optical and IR emission lines 
throughout the orbital period and find evidence, in the Balmer lines, for an optically 
thin accretion disk and an optically thick hot-spot. Optical and IR emission lines are 
used to compute radial velocity curves. Fits to our radial velocity measurements give 
an internally inconsistent set of values for Kl, 7, and the phase of red-to-blue cross- 
ing. We present a probable explanation for these discrepancies and provide evidence 
for similar behaviour in other short orbital period dwarf-novae. Selected optical and 
IR spectra are measured to determine the accretion disk radii. Values for the disk 
radii are found to be strongly dependent on the assumed WD mass and binary orbital 
inclination. However, the separation of the peaks in the optical emission line (i.e. an 
indication of the outer disk radius) has been found to be constant during all phases 
of the super-cycle period over the last 40 years. 
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1 INTRODUCTION 

Dwarf-novae, among the cataclysmic variables (CVs), are 
semi-detached binary star systems where a white-dwarf (pri- 
mary star) accretes gas from a cool companion (secondary 
star). The gas from the secondary star accretes through the 
inner Lagrangian point (LI), into the primary Roche lobe 
and slo wly spirals in toward the white-dwarf. The in-flowing 



emitted flux distribution over the disk, determines the line 
profile shape. Consequently, spectral line analysis is a funda- 
mental method employed to understand the physical state 
of the accretion disk and its spatial flux distribution. The 
most common emission lines used in such analyses are those 
due to hydrogen and helium. 

In this paper, we continue the analysis of the accre- 



tion disk in WZ Sge which was started in skidmore et al 



gas forms an accretion disk which orbits around the white- 



dwarf before actually being accreted onto the primary star. 
During the time the gas is within the accretion disk, current 
theory tells us that half of the gravitational energy of the gas 
is radiated outward while that which remains heats the ac- 
cretion disk. Indeed most of the visible and IR light output 



200C (hereafter PAPER I). PAPER I presented Doppler to- 



by a typical dwarf-novae is from the accretion disk (Ciardi 



et al. 1998, Warner 1995). Emission lines seen in the spec- 



tra of a dwarf-nova arise from the vicinity of the accretion 
disk and can be double peaked due to the Doppler m otions 
of the disk within the binary ( Home fc Marsh 198(j ). The 



velocity of the material within the disk, together with the 



mography of optical and IR disk emission lines as well as 
ratioed Doppler maps. General conclusions reached there in- 
dicated different optical depth at the hot-spot and the accre- 
tion disk, and an accretion disk structure dissimilar to those 
predicted by models. In this paper (PAPER II), we perform 
spectral analysis by direct measurements of the emission line 
profiles. We organize the present paper in three main sec- 
tions. Section 2 describes the orbital phase dependence of 
Ha, H/3, and IR emission lines. Section 3 presents newly de- 
termined radial velocity curves using optical and IR emission 
lines. We show that the amplitude and phase zero of the ra- 
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Table 1. Log of the observations of WZ Sge used in this analysis. Velocity dispersion is stated for the line used; it is not the mean 
velocity dispersion of the complete data set. (t average value for all three lines). 



Date of 
observations 


Wavelength 
range (A) 


Emission 
Lines 


Dispersion 
(A/pixel) 


App. dispersion 
(km sec -1 ) 


Exposure time 
(sec) 


Run duration 
(hours) 


Orbits 
covered 


Number of 
spectra 


27/07/96 


6375 - 6780 


Ha 


0.40 


18.1 


40 


5.51 


4.1 


237 


27/07/96 


4590 - 4995 


Hj9 


0.40 


24.6 


40 


5.03 


3.7 


239 


28/07/96 


6375 - 6780 


Ha 


0.40 


18.1 


40 


4.07 


3.0 


139 


28/07/96 


4590 - 4995 


H/3 


0.40 


24.6 


40 


4.10 


3.0 


137 


28/05/97 


18000 - 24600 


Br7, Br5, He I 


13.2 


193f 


480 


2.70 


2.0 


18 


08/08/98 


10300 - 13400 


Pa/3 


6.22 


145.6 


240 


6.73 


3.5 


30 
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Sample Ha and H/3 phase binned spectra 
-2 „-l 



(50 bins per orbit). Each spectrum is offset from the one below by 2.1 



Figure 

10 -15 erg cm -2 s — 1 A -1 and each bin is labeled with mid-orbital bin phase based on the optical photometric eclipse ephemeris listed in 
PAPER I. Note that the Ha spectra show a weak He I emission line (A 6678A), which follows the Ha S-wave motion (see trailed-spectrum 
in figure 1 of PAPER I). The narrow "absorption features" seen in the H/3 spectra are due to bad pixels. Note: the left and right panels 
have the same flux units. 



dial velocity curves are systematically biased. We show that 
this is due to different gas opacities in the hot-spot and the 
accretion disk, the hot-spot having a larger opacity than 
the accretion disk. We further develop the idea of increas- 
ing hot-spot contamination in systems where the difference 
between the optical depth in the hot-spot and the accretion 
disk is large, these systems being short orbital period, low 
mass transfer rate dwarf-novae. In section 4 we present our 
determinations for the accretion disk radii and compare our 
values with previously reported determinations. 



phase bins per orbit in order to improve the signal-to-noise 
ratio of the data. 



2 EMISSION LINES 

We analyze the emission lines evolution by inspecting their 
profile, measuring the fluxes and computing the Balmer 
decrement (last analysis on optical spectra only). Figure ^ 
summarizes the results from each analysis. 



Complete details concerning the data set used here (see 
table |l|), are in PAPER I. All data are phased according 
to the photometric eclipse ephemeris stated in table 7, PA- 
PER I. Optical spectra of WZ Sge were grouped into 50 



2.1 The line profiles 

Figures |l| and [| present a sample of optical and IR spectra 
respectively to show the complex emission line profiles and 
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2.1x10 4 2.2x10 4 
Wavelength (A) 




1 .2x1 4 
Wavelength (A) 



Figure 2. IR spectra in the K and J bands. Spectra are labeled with the phase at mid-exposure time as set by the photometric eclipse 
ephemeris listed in PAPER I. Flux is in mjy and the offset between each spectra is 2.8 and 4.0 mjy for K and J spectra respectively. Note: 
the left and right panels have the same flux units. 



their dramatic evolution throughout an orbital period. Both 
optical and IR emission line profiles evolve similarly and 
owe most of their changes to the variable strength of the 
hot-spot around the orbit. Phases 0.3, 0.5, 0.84 and 0.63 are 
particularly suitable to point out the hot-spot evolution over 
the orbit. 

Inspection shows that phase 0.3 is when the hot-spot, 
in both Ha and H/? emission lines, is at the center of the 
accretion disk emission, i.e. has zero radial velocity with re- 
spect to an observer co-moving with the white-dwarf. From 
the system geometry determined in PAPER I we expect the 
hot-spot in the middle of the accretion disk emission line 
at phase ~0.3 only if the gas at the stream-disk impact re- 
gion has just the stream velocity (see figure ^). This con- 
dition implies that the accretion disk density is low enough 
such that the gas at the stream-disk impact region keeps its 
stream velocity with little effect from the Keplerian motion 
of the underlying accretion disk material. In the IR spectra 
the phase of the hot-spot zero radial velocity occurs slightly 
later, i.e. at phase 0.33. Phase binning of the optical spectra 
with the same time resolution of the IR spectra has showed 
that such a delay is significant and implies different velocity 
of the hot-spot gas in different wave-bands. We explain such 
a phase delay with the hot-spot emission from different lines 
spreading along the stream trajectory (figure ^). 

The hot-spot appears again at the center of the emission 
line at phase 0.84, both in the optical and IR spectra. The 
expected phase for the optical spectra is 0.80 in the assump- 
tion of symmetric and isotropic hot-spot emission. Further 



evidence for asymmetric and anisotropic hot-spot emission 
is in the hot-spot profile which appears double peaked at 
phase 0.84. This feature is real as it is common to both Ha 
and H/3 binned spectra and is visible for a phase interval of 
<0.5 (see also the trailed spectrograms in figure 1 of PA- 
PER I). A detailed study of the double structure is beyond 
the scope of this paper. The lower spectral resolution in the 
J and K bands does not resolve such a feature if present (the 
hot-spot peak separation is ~9A in the Balmer lines). 

Phase ~0.53 is when we expect the maximum hot-spot 
blue-shift and enhanced blue peak emission in the accretion 
disk lines. Observations do not confirm this expectation and 
show red and blue shifted peaks of roughly equal strength 
between phase 0.5 and 0.56, both in the optical and IR spec- 
tra. We explain this observation with little or no contribu- 
tion from the hot-spot to the overall emission line flux. This 
is also evidence of the anisotropic hot-spot emission. 

Asymmetries in the accretion disk are evident in the 
Ha emission line which displays a stronger red-shifted peak 
at phase 0.63. Doppler maps in PAPER I show asymmetric 
accretion disks also in the H/3, Br7, He I, and Br<5 lines. The 
accretion disk emission line profile evolves from a shallow 
U-shaped profile to a V-shaped profile, between phase 0.5 
and 0.63. This is particularly evident in the optical spectra; 
among IR observations, only the J band spectra show similar 
evolution. Home fc Marsh (1986 ) predict U-shaped emission 
line profiles for optically thin accretion disks, and V-shaped 
profiles for the optically thick case. This might suggest that 
the changing profile of the emission lines is evidence that 
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EMISSION LINE 



CONTINUUM 



H/S 



He 



ivy 



- 1.2x10 



Phase Phase 
Figure 3. Left: the line flux for each emission line vs the orbital nhase; Right: the continuum flux underlying the emission lines vs orbital 
phase. Fluxes measured at phase 0.3, 0.5, and 0.84 are in table El Both emission line and continuum fluxes were measured using the e-e 
command in IRAF task splot; fitting was not attempted. Average uncertainties on measurements are plot on the first data point in each 
panel; error-bars on Ho, Br7, and He I 2.06/xm emission lines are smaller than the data point size. Note: the left and right panels have the 
same flux units. 



Line 



Max Flux 



Min Flux Flux at 4> ~ 0.3 



Ha 
H/3 
Br7 
He I 2.06/mi 



11.00xl0~ 14 
4.95xl0~ 14 
1.19X10" 14 
5.67xl0~ 15 



6.75X10" 14 
2.65X10" 14 
l.OlxlO -14 
3.54xl0~ 15 



8.80X10" 14 
4.00X10" 14 
8.32xl0~ 15 



Table 2. The maximum and minimum fluxes, and the flux near 
phase 0.3, as measured in each emission line. Flux units are erg 
cm 2 sec -1 A -1 . The minimum fluxes are observed about phase 
0.5 and are believed to correspond to the phase when there is no 
contribution from the hot-spot. The maximum flux is observed at 
phase 0.84 when the hot-spot is viewed from the outside. Phase 
0.3 is the phase when the hot-spot is viewed from the inside. 



some areas of the disk have optically thin line emission and 
some have optically thick. 



2.2 Line flux 

The emission line flux was measured by manually marking 
the position of line wings with the cursor and summing up 
the flux from each pixel in between; the continuum flux cor- 
responding to a linear interpolation of the same two cursor 
position was automatically subtracted during the line flux 
measurements. Each of the emission lines measured (with 
the exception of the Pa/3 line) show quite strong modulations 



apparently correlated with the line profile evolution (fig- 
ure^. The modulations display maxima at phase 0.84 and 
0.3, and a minimum at phase 0.5. Thus, the hot-spot shows 
stronger emission when viewed from the outside (phase 0.84) 
and from the inside (0.3), while it does not contribute signif- 
icantly to the emission line flux when viewed from upstream 
(phase 0.5). Assuming that the flux at phase 0.5 is due only 
to the accretion disk (see figure |l] and figure ^, and previous 
subsection) we can compute the relative contribution by the 
hot-spot to the total emission line flux. At phase 0.84 the 
hot-spot emits 50% of the total line flux in H/3 and ~ 30% 
in Ha and the IR lines. When, viewed form the inside the 
hot-spot relative contribution decreases to 38% in H/3 and 
to ~ 20% in Ha and Br7. The He I 2.06/im line does not dis- 
play an evident hump at phase 0.3, thus the relative hot-spot 
contribution was not computed there. Assuming a spherical 
hot-spot with radius S = 3.1 x 10 cm (fcmak 19931 ), cylindri- 
cal accretion disk with thickness 2H ~ 3.8 x 10 cm and 
inclination i — 75° (Smak 1993), the decrease of the hot-spot 



flux from phase 0.84 to phase 0.30 due to occultation by an 
optically thick disk is 70%. We observe a decrease of only 



The value was derived comp uting the scale height H as given 
in ( Williams fc Shipman 198q ) assuming a uniform temperature 
T=6000K in the accretion disk and a white dwarf mass of Mvkd = 
0.82M Q (PAPER I). 
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to a white-dwarf mass of 0,82Mg ( PAPER. I), an d the secondary 
star mass (0.058Mq) computed by rimak (199s| ). Units on both 
the x and y axis are m. Straight full lines mark the line of sight 
at the phases considered in the text. Dashed and dotted lines 
show the directions of stream and Keplerian velocity respectively 
for the illustrative case of Ho hot-spot position. Crosses along 
the stream trajectory qualitatively mark the hot-spot position for 
each emission line as predicted by our results. From the left to 
the right: Ha, H/3, Pa/3, and Br7 in order of increasing excitation 
energy and gas temperature. 



30% in H/3 and ~50% in Ha and Br7. Thus either the accre- 
tion disk gas is at least partially transparent to the hot-spot 
emission line flux, or the assumed geometry inappropriate. 

The continuum light curve is expected to be double- 
humped as already observed by Ciardi et al. (1998) in the 



IR, and by Patterson et al. (199' 
the optical. Right panel of figure 



and 



Skidmore (1990 in 



shows that only our K- 
band continuum measurements match with previous obser- 
vations. The optical continuum shows large uncertainties on 
each data point and poor modulation throughout the or- 
bit as the method used to measure line and continuum flux 
does not take into account the underlying white-dwarf ab- 
sorption. The IR spectra, which do not show any evidence 
of white-dwarf features have the expected double-humped 
structure and smaller uncertainties in the continuum light 
curves. Fluxes from the Pa/3 emission line and its underlying 
continuum have not been reported in figure ^ and table g 
because the data show variations throughout the run of ob- 
servation. The run covered about four cycles and variations 
between cycles do not appear to be correlated. We cannot 
say whether this is due to intrinsic variability of the emitting 
source or to external causes. 



Figure 5. The Balmer decrement, D v (Hcti 'Hf3), versus orbital 
phase both for the accretion disk (top) and hot-spot (bottom). 



2.3 The Balmer decrement 

We measured the Balmer decrement for both the accretion 
disk and the hot-spot region to check for different opacities. 
The Balmer decrement is defined by the ratio of the line in- 
tensity in frequency units D v (Hot/ Hf3) = I v (Ha) / I V (H f3). 
The Balmer decrement, as measured from spectra, relies on 
line fluxes, because the line intensity is broadened by effects 
due to both the intrinsic properties of the emitting gas and 
the instrumentation used during the observation. The com- 
parison of measured flux-ratios to tabulated intensity-ratios 
assumes that the emission lines of interest have identical pro- 
files, each broadened in a similar way. We observe Balmer 
lines which are far from being identical in shape, however, 
table 6 of PAPER I show that time averaged disk profiles 
have identical velocity broadening. De-blending of the com- 
ponents is not possible without a complex model of hot-spot 
shape and accretion disk structure, thus we measured the 
peak intensity for each component, on continuum-subtracted 
spectra. We measured the disk emission at each phase using 
only the red or blue line peak, whichever was unaffected by 
the hot-spot emission at the given phase. The hot-spot line 
intensity has always some underlying disk emission blending, 
except near phase ~ 0.3 and ~ 0.84, when the hot-spot emis- 
sion is at the center of the emission line and the relative con- 
tribution from the accretion disk is negligible. The top panel 
of figure |E| shows that the Balmer decrement within the ac- 
cretion disk is approximately constant throughout the orbit 
with an average value of D v (Ha/H/3) = 3.82, which is larger 
than the predicted v alue of D v (Hal H 0) ~ 1 for an optically 



thick accretion disk ( [Williams 198C| , [Williams 1983[ ). I v (Ha) 
and I v (H0) each show phase dependent modulations, im- 
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plying a non-uniform accretion disk emission, but their ra- 
tio remains constant. However, despite the uncertainties, 
the Balmer decrement in the accretion disk shows a clear 
trend around phase 0.3 toward higher opacities. The bot- 
tom panel of figure [B| shows the hot-spot Balmer decrement 
vs orbital period. The Balmer decrement here is strongly 
phase dependent and shows two deep minima at phase inter- 
vals 0.2-0.3 and 0.8-0.84, i.e. when observing perpendicular 
to the stream trajectory at the hot-spot position and mea- 
suring the hot-spot line intensities without any underlying 
accretion disk emission. The hot-spot Balmer decrement is 
D v (Ha/H(3) ~2. The value of D„(Ha/H/3) outside phases 
0.2-0.3 and 0.8-0.84 is ~ 3 and approaches the average ac- 
cretion disk Balmer decrement when the hot-spot emission 
is weak. In the hot-spot, I„(H/3) shows a stronger modula- 
tion around the orbit than I v (Ha) which is largely biased 
by the underlying accretion disk emission. 

The Balmer decrement provides information about the 
gas temperature and density where the emission line forms; 
it gives a unique gas density for each assumed temperature. 
However, care must be taken because different temperatures 
may give the same value of the Balmer decrement depending 
on the assumption of optically thick or optically thin emis- 
sion lines. In the interpreta tion of our results we mainly use 
models by Williams (1991) who modeled emission lines in 



optically thin accretion disks. He computed H/3 strengths 
and Balmer decrements for a grid of temperatures, orbital 
inclinations, and mid-plane nucle on densities. In figure^ we 
plot the results of Williams (1991) together with the Balmer 
decrement we measure at the hot-spot (lower horizontal line) 
and the accretion disk (upper horizontal line) in WZ Sge. 

The value of D v (Ha/ H /3)=2 for the hot-spot region 
corresponds to a black-body optically thick emission with 
a temperature of T<5000K. This decrement apparently 
matches a number of observed Balmer decrements in CVs 



(Williams 1983). However, because we observe the He I emis- 
sion line at 2.06^im, parts of the hot-spot region must be at 
least 15000K (see section An optically thin gas with a 
tempertiture of 15000K and a Balrher decrement of 2 cotte- 



sponds Lu a gas density uf lug N — 11.0 in Williams' mudel 

(see figure |), which is lower than the derived gas density 
within the accretion disk (see below). The gas density in 
the stream is expected to be 3 to 4 t imes larger than th e 
gas density within the accretion disk ( Lubow fc Shu 1976 ). 
We thus conclude that the hot-spot is at least partially op- 
tically thick and characterized by a steep temperature gra- 
dient. The Balmer emission lines form down-stream from 
the initial impact region from a relatively low temperature 
(T~5000K) optically thick gas, and the He I 2.06pm emis- 
sion line forms in the outer stream-impact region where the 
temperature is higher (see section |^) . 

The Balmer decrement D v (Ha/H/3) = 3.82 of the gas 
within the accretion disk would correspond to a black-body 
of temperature T~3400K. However, we expect the temper- 
ature to be higher as in an optically thick a-disk the tem- 
perature is never expected t o fall below 600 0K over a wide 
range of mass accretion rate (Williams 1980). Moreover, Ra- 



tioed Doppler tomograms and radial disk profile of WZ Sge 
(see PAPER I) suggest a partially optically thin accretion 
disk. We conclude that the accretion disk is optically thin 
in the Balmer emission lines. According to Williams (1991) 
a Balmer decrement of D v {Hct/ H (3) = 3.82 corresponds 
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Figure 6. Ba lmer decrement, D v {Ha/ H (3), as predicted by 
Williams (1991) for optically thin disk emission lines. Each line 



represents the Balmer decrement as function of the mid-plane 
accretion disk density, N cm -3 , for a given temperature and or- 
bital inclination. The models for T=8000, 10000, and 15000K, all 
approach D v (Ha/ H 0) ~ 1 at high density (log N ~ 13. — 13.5), 
as expected for an optically thick gas. The horizontal lines 
mark the Balmer decrement measured for the accretion disk 
{D v {Ha/Hf3) = 3.82) and the hot-spot (D v {Hct/Hf3) = 2) in 
WZ Sge. 



to an optically thin gas of temperature 6000K and den- 
sity logiV ~ 13.5 (see figure |^). A roughly similar value 
of log N e ~ 1 3 for gas at a temperat ure of T =5000K may be 



derived from 



Drake fe Ulrich (1980|) models. [Drake & Ulrica 



(1980| ), using the escape probability approach, determined 



Balmer decrements and hydrogen lines strength for an infi- 
nite slab of gas at various opacities, temperatures, and den- 
sities. We can use these estimates to determine the optically 
thin accretion disk mass assuming: i) a uniform gas density 
of log N — 13.2, ii) a cylindrical accretion disk with radii as 
determined in PAPER I (see also table ^), and height 2H as 



already defined in subsection 2.2. The result is an accretion 
disk mass of Md ~ 5.7 x lO _lb M0, which is about 5 orders 
of magnitude smaller than the total ma ss accreted during 
the inter-outburst period (~ 1O~ 9 M0, [5mak 1993 ). It is 
also ~4 order of magnitude smaller than the critical mass 
(i.e. the maximum accretion disk mass before an outburst) 
expected in low mass transfer rate systems However, typ- 
ical values for the mass transfer rates and critical masses 
have been computed using the assumption of a-disk models 



t We derived the value M 



crit ~ 4x 10 11 Mq either using average 



val 



les for the accretion disk ra,dii a 



d surface density as predicted 
by |Lasota, Hameury fc Hure (1995] ) (Ri n ~ 4.5 X 10 9 cm flout ^ 



Q lu cm. logS ~ 2.4, for M ~ 10~ n M^yr" 1 ), or formula 
f!8l inpmak Q993|) with a ~ 0.0 03 as predicted for TOADs by 
Howell,*Szkody fc Cannizzo (1995| ). 
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which are optically thick and cannot t herefore be comp ared 
directly with optically thin gas models. Williams (1991) pre- 



dicts an accretion disk that is only partially optically thin 
at a mass transfer rate of M ~ lO _11 M0yr _1 , the optically 
thin region starting from a radius of 5 x 10 9 cm, outward. 
This value seems in quite good agreement with the radius 
of r ~ 2.6 x 10 9 cm at which we observe the transition of the 
Balmer line ratio from optically thick to optically thin case 
(see section 4 and figure 7 in PAPER I). 



3 RADIAL VELOCITY CURVES 

Optical and IR spectra were used to measure radial velocity 
curves for H/3, Ha, Pa/3, Br7 and He I (A2.06/mi) emission 
lines. Radial velocity curve measurements were not produced 
for the two emission lines He I Al.OS^m and Pa7 as they were 
strongly blended together. 

Figure [7] presents our radial velocity curves and their 
best fit for the five emission lines listed above. Each data 
point corresponds to the radial velocity of the emission line 
central wavelength. Central wavelengths were computed by 
using a technique similar to the Pogson's method to deter- 
mine the m id-time of an eclipsing binary with an asymmet- 
ric eclipse (Glasby 1970). The radial velocity measurements 
were fit by a grid-search method minimizing x 2 ■ The sinu- 
soidal fitting function used was: 



■y + KiX sin[27r((/>i + 0.50 - <f> a )] 



(1) 



where Vi is the measured velocity and <f>i the observed pho- 
tometric phase. The fitting parameters 7, Kl, and (f> are 
respectively the systemic velocity, the Keplerian velocity of 
the white-dwarf around the system's center of mass, and the 
red-to-blue crossing in the radial velocity curve, correspond- 
ing to secondary inferior conjunction. We added the phase 
shift of 0.50 in Eq.|l] in order to have <f> exactly matching 
the secondary star inferior conjunction, and not the primary 
star inferior conjunction. 

Our best fit parameters from each emission line repre- 
sent an inconsistent set of values (see table ^|) . The systemic 
velocity, 7, runs from negative to positive values, spanning 
a range of ~ 300km sec -1 wide. The white-dwarf Keple- 
rian velocity, Kl, runs from ~46 to 120 km sec -1 ; while the 
red-to-blue crossing, <f> , swings about phase 0.19, covering 
a range A(j> ~0.1 in phase. Interestingly, a similarly incon- 
sistent set of system parameters was determined for the op- 
tical radi al velocity curves of the SU UMa type dwarf-nova 
VY Aqr (|Augustcijn 1994|). 



When fitting and interpreting the sinusoidal radial ve- 
locity curves, it is assumed that: i) the disk orbits with Kep- 
lerian velocity around the white-dwarf, and ii) the accretion 
disk follows the white-dwarf motion around the binary's cen- 
ter of mass. These approximations fail when hot-spot emis - 
sion and disk asymmetries are not negligible, smak (1970), 



shows that the hot-spot can induce a measured Kl larger 
than the true value and a spuriously delayed moment fo r 
inferior conjunction. Paczynski, Piotrowski & Turski (1965) 



explain a larger Kl amplitude with non Keplerian motion. 

We believe that our measured radial velocity curves re- 
sult from combination of two different motions: 1) the mo- 
tion of the accretion disk as described at points i) and ii) of 
the previous paragraph; and 2) the motion of the gas within 




12.4 12.6 12.8 13 

Excitation Potential Energy (eV) 

Figure 8. The linear relationship between the observed phase off- 
set and the excitation potential of the emission line. The straight 
line is just to highlight the linearity of the observed distribution. 
The slope is 0.088(phase/eV); the y-intercept is -0.93(eV). 



the hot-spot region which moves along the ballistic trajec- 
tory at the stream velocity. Combination of the two velocity 
components in points 1) and 2) will change the amplitude 
and phase of a sinusoidal radial velocity curve, the largest 
departures occurring in systems where the hot-spot emission 
represents the largest fraction of the total emission line flux. 
An anisotropically emitting hot-spot may affect the 7 veloc- 
ity measurements. Optically thin accretion disks are likely 
to be low density gas regions (see section 3.1), which neither 
affects the stream velocity of the in-falling gas, nor absorbs 
the flux emitted at the impact region. Thus, emission line 
profiles from optically thin accretion disks are expected to 
be heavily biased by the hot-spot emission. 

The phase offset is found to be correlated with the en- 
ergy required to produce each emission line. Figure ^ shows 
the observed phase offset versus the excitation potential of 
each H emission line as given by the standard formula. Fig- 
ure ^ clearly shows a linear relationship between the ob- 
served phase offset and the emission line energy and implies 
a smeared out hot-spot region with decreasing excitation en- 
ergy and temperatures as the gas moves down-stream. Re- 
sults si milar to those in figure [S| may be obtained using the 
data in 



Augusteijn (1994) on VY Aqr. We fit the phase offset 



computed by Augusteijn for three hydrogen emission lines 
and found a linear relationship with almost identical slope 
and y-intercept as in figure |§|. 

The He I 2.06^im emission line was not considered in fig- 
ure ^. The He I line 2.06/im is a transition strongly coupled 
with the resonant line at 584 A. Both lines decay from the 
same 2 X P level but the 584A transition has a probability of 



occurrence which is ~ 10 times larger (Najarro et al. 1994) 



Strong He I emission at 2.06/xm is rarely seen in astronom- 
ical objects and observed only in the extended atmosphere, 
wind and/or disk- like formations of WN and Be super-giant 
stars, where they form because of the high temperature, 
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Figure 7. The radial velocity measurements and their best fit sine curves. Phasing is based on the photometric eclipse ephemeris of 
Skidmore (199£), the period was fixed at the orbital period, while the amplitude, Kl, and the phase <j> were free parameters in the 



fit. The red-to-blue crossing of each line is shown as a vertical dashed line. The mean wavelength or systemic velocity, 7, is shown as a 
horizontal dotted line. 



Table 3. Results of fitting the radial velocity curves (See Fig [?]). 



Line Rest Wavelength V/R crossing K\ amplitude Systemic velocity 
A (j> km sec -1 (7) km sec -1 



H/3 


486L33 


0.1910±0.0058 


82.0±2.5 


-84.6±2.0 


Ha 


6562.7725 


0.1331±0.0061 


67.9±2.8 


-71.4±1.9 


Pa/3 


12818.1 


0.246±0.038 


46.2±8.5 


10.4±7.0 


He I 


20581.3 


0.157±0.014 


121.4±12.6 


144.1±9.2 


Br7 


21655.3 


0.228±0.023 


108±15 


208±11 
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Star Name 



DN Type Orbital Period Phase Offset 



Kef. 









L -* ( r — ro <rphot 


(dee) 




OY Car 


SU 


.06312 


.1167± .0058 


83.3±0.2t 


(1) 


HT Cas 


SU 


.073647 


.095± .009 


76.4 


(2) 


Z Cha 


SU 


.0745 


.0904 ± .0017 


81.7 


(3) 


EM Cyg 


zc 


.29090950 


.057 ± .011 


63±10t 


(4) 


U Gem 


UG 


.17708 


.015 ± .011 


69.7±0.7t 


(5) 


V2051 Oph 


UG? 


.062427887 


.132 ± .020 


80.5+2+ 


(6) 


CN Ori 


UG 


.163190 


.03 ± .01 


65^70 


(7) 


IP Peg 


UG 


.15820 


.067 ± .008 


79.3±0.9 


(8) 


WZ Sge 


SU 


.05668684707 


.1620 ± .0042 


75±2 


(9) 


DV UMa 


SU? 


.0859722 


.10 ± .02 


72 


(10) 


HS 1804+6753 


UG 


.2099370 


.044 ± .002 


84.2±0.6 


(11) 



Tabl e 4 . Orbital period, phase offset , and inclination for the sample of dw arf-n ovae, see figure^. References: (1) jjchocmbs fc ilartma~ 
(1983^1 ; |Bailey fc Ward (198l|l: (2) Young. Schneider fc Shcctman (108l|'l: [3] [Marsh. Home" Shiuman (1987] ); (4) |Stovcr, Robinson 
& Nath cr (1981 )77^ Smak (1976); (6) Watts ct al. (198Q; (7) |Barrcra fc Vogt (198^ ) ; (8) Vlarsh (1988 ); (9) this paper; (10) Szkody & 
Howell (1993[): (11) |Ficdlcr, Barwig fc Mantel (1997|). t fitter fc Kolb (1998|)" 



velocity and density of the out-flowing gas (15000- 30000K, 
~700 Km s ec~ , and ~ 10 _J Mq yr _1 , respectively, Najarro 
et al. 1394). The presence of the He I 2.06/im emission re- 
quires the He I 584A line to form in an optically thick, high 
temperature 15000-30000K region, which, in short period 
dwarf-novae, can only occur in the shock heated stream-disk 
impact region. Formation of He I 2.06/^m in a hot corona is 
excluded by the Doppler tomogram in PAPER I figure 3. 

We interpret our results above to provide a picture for 
the quiescent accretion disk in WZ Sge (figure ^) , where each 
hot-spot emission line is likely to be aligned along the stream 
trajectory into the accretion disk, with the higher energy 
emissions arising at the higher temperature outer edge of 
the accretion disk impact region, and the lower energy lines 
primarily forming further down-stream. 



3.1 Phase offset as a measure of the hot-spot bias 

White-dwarf and/or hot-spot eclipses are present in the 
broad band photometric light curves of many dwarf-novae. 
In the case of the white-dwarf eclipse, the observed eclipse 
minimum defines the time of the secondary star inferior con- 
junction. In systems displaying only a bright spot eclipse, 
the phase offset between binary phase zero and the hot-spot 
eclipse is small (<0.05 in phase, Hantzios 1988), so that we 
still have a good idea of the time of the true binary phase 
zero. Similarly, when measuring radial velocity curves, we 
expect the spectroscopic ephemeris to match the photomet- 
ric one in systems with perfectly symmetric accretion disk, 
and a phase offset larger than zero in systems having asym- 
metries in their disk. In particular, a larger phase offset is 
expected in systems where the hot-spot represents the main 
asymmetry in the accretion disk and a considerable fraction 
of the total emission line flux. This is the case of WZ Sge 
given that our observations average out possible short term 
variation and Doppler tomograms in PAPER I show accre- 
tion disk asymmetries always considerably fainter than the 
hot-spot emission; moreover in section ^| we showed the ac- 
cretion disk to be optically thin. To test the hypothesis of 
phase offset as a measure of the hot-spot bias and the op- 
tical properties of the accretion disk gas we searched the 
literature for measured dwarf-novae phase offsets and plot 



, DN, SU Type 

x DN, UG Type 

t DN, ZC Type 



0.05 0,1 0.15 0,2 0.25 0, 

Orbital Period (days) 

Figure 9. The observed phase offset vs. orbital period. Phase off- 
set is the delay of the spectroscopic secondary conjunction (red 
to blue crossing) as measured from the disk emission lines, with 
respect to the photometric eclipse minimum (i.e. the photomet- 
ric phase zero). Different dwarf- novae subtypes are plotted with 
different symbols. Dotted vertical lines delimit the period-gap. 



these values against the binary orbital period. Theory pre- 
dicts low mass transfer rates in short orbital period systems 
(Howell, Rappaport & Politano 1997), implying low density 



possibly optically thin accretion disks. 

Phase offsets have been observed in 



several different 



CVs (Watts et al. 1986 



Shaffer 1991 



and Stover 1981), but 



until now th ere has been no ex plicit or systematic search 
for its cause. Watts et al. (1986), observed phase offsets in 
both nova-like and SU UMa stars, claiming that the hot - 
spot cannot be the cause of such an offset. [Stover (1981 ) 



concluded that the best candidate to explain the phase off- 
set observed i n five dwarf-novae is only the hot-spot. Finally, 
Shaffer (1991 ) collected most of the observed phase offsets in 
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the literature and plotted them vs. the binary orbital period 
but did not provide a possible explanation. 

The eleven dwarf-novae listed in table ^ have been 
carefully selected in order to provide a homogeneous sam- 
ple. This was achieved by selecting: i) high inclination 
dwarf-novae that show photometric eclipses, ii) only Balmer 
line radial velocity curves, and Hi) publications reporting 
both photometric and spectroscopic analysis. Restrictions 
for point i) are to avoid biases other than the hot-spot and 
to reduce unknown inclination effects. Point ii) reduces the 
scatter in the measured phase offset by focusing on emis- 
sion lines arising from regions with similar physical condi- 
tions, and finally, point Hi) guarantees consistently measured 
phase offset and reduces any bias due to possible changes 
in the quiescent accretion state. Figure [] plots our selected 
sample which provides evidence for the expected correlation: 
the phase offset is larger in short orbital period systems. 



4 ACCRETION DISK SIZE 

It is possible to estimate the inner and outer accretion disk 
radii in a dwarf-nova by measuring the emission line wing 
and peak separations respectively, and assuming a Keplerian 
velocity. For an annulus of gas in Keplerian motion at a 
distance r from a white-dwarf of mass Mwd , orbiting at an 
angle i to the line of sight, the Keplerian velocity is: 



Half Wing Separation (Km sec 1 ) 



(2) 



In a double peaked emission line from an accretion disk, 
the Keplerian velocity of the outer disk radius corresponds 
to half the line peak separation (hereafter HPS); while, the 
Keplerian velocity at the inner accretion disk radius corre- 
sponds to half the line wing separation (hereafter HWS). 

It is important both to have reasonable system parame- 
ters and to measure HWS and HPS on unbiased spectra, i.e. 
spectra in which the wings and peaks of the accretion disk 
emission line are not affected by hot-spot emission (figure |l| 
phase 0.3 and 0.84). Spectra in which the hot-spot emission 
clearly contaminates wings and/or peaks of the accretion 
disk emission line profile may provide a measure of the hot- 
spot bias. Both HWS and HPS were measured reading the 
cursor position with the IRAF task splot. HWS values from 
each spectra are the average over 5 measurements. 

4.1 Measurements on uncontaminated spectra 

We selected a sample of uncontaminated phase binned Ha 
and H/3 spectra. The assumption is that the hot-spot emis- 
sion does not affect the wings of the accretion disk emis- 
sion lines. This selection criteria produced three binned Ha 
spectra (at phases 0.28, 0.30 and 0.32) and six binned H/3 
spectra (at phases 0.28, 0.30, 0.32, 0.34, 0.82 and 0.84). We 
applied the same selection criteria to IR spectra and selected 
two uncontaminated spectra in both the J and K band (at 
phases 0.33 and ~0.84 in each case). In case the hot-spot 
feature was not evident (figure ^, J-band, phase 0.3) we se- 
lect those spectra which are the closest in phase to the un- 
contaminated optical spectra listed above. HWS and HPS 
were first measured on each bin separately, and then aver- 
aged. HWS and HPS measured values are in table ^[ We see 



Uncontaminated Spectra? 



Line 



Averaged Spectra 
Kl subtracted 



Ha 
H/3 
Pa7 
He I (2.06/im) 
Br7 



1606 ± 65 
1692 ± 98 
1868 ± 260 
1740 ± 97 
1598 ± 50 



2065 ± 7 
1663 ± 23 
2063 ± 75 
1999 ± 75 
1601 ± 20 



Table 7. Half wing separation value corrected for the white-dwarf 
orbital motion. Each emission line measure has been corrected 
with the Kl velocity determined in the corresponding radial ve- 
locity curve, t Same values as given in table 5l 



that values from optical and IR emission lines are in good 
agreement each other, the exceptions being the He I 2.06/^m 
and Pa/3 lines which apparently have smaller HPS value and 
larger HWS values than in the Ha, H/3, and Br7 measure- 
ments. We applied a t-student test (to a confidence level of 
a = 0.05), to check for statistical differences between pairs 
of values. The test showed that the HWS values have to be 
considered statistically identical in all the emission lines, i.e. 
the same inner disk radii; while the H/3 HPS is statistically 
larger than the measured HPS from Ha, Br7, He I 2.06/im, 
and Pa/3, thus implying a smaller outer disk radius (~ 10%). 

Velocities from uncontaminated spectra were used in 
equation ^ to compute the accretion disk radii for a variety 
of primary masses and inclinations (see table Q). It is evident 
that small white-dwarf masses (Mwd < O.5M0), imply a 
small accretion disk outer radius (< 50% of the primary 
Roche lobe radius) and a relatively large ratio Rwn/Rm 
(< 0.5); while larger masses (Mwd > 1M©) give a large 
outer disk radius (up to 90% of the primary Roche lobe 
radius) and Rwo/Rm ~ 0.07, implying an annulus shaped 
accretion disk. The input value of the white-dwarf mass does 
not affect the ratio, R=Ri„/R ou t, which is equal to 0.2 in 
each case. 



4.2 Measurements on the averaged spectra 

Optical and IR average spectra show relatively symmetric 
accretion disk emission lines. This is expected because the 
hot-spot motion, i.e. the S-wave, smoothes out the hot-spot 
features at each phase, and evenly distributes its emission 
throughout the orbit. Average spectra also show broader 
emission lines than single spectra. This happens for two rea- 
sons: i) the common orbital motion of the accretion disk 
and the white-dwarf around the binary center of mass will 
broaden the disk emission line by a term Kl, i.e. the white- 
dwarf Keplerian velocity; and ii) the wings of the hot-spot 
emission blend with the blue and red-shifted accretion disk 
emission lines. The S-wave velocity is usually smaller than 



the outer disk edge Keplerian velocity (Gilliland, Kemper 



Suntzeff 1986; hereafter GKS), thus, we expect average spec- 



tra with broader accretion disk emission peaks and smaller 
peak separation than observed in single spectra. HPS and 
HWS from average spectra are given in table |^. HPS val- 
ues from optical and IR emission lines agree fairly well and 
confirm -with the exception of He I 2.06/im and Pa/3- the 
expectation. HWS values are in good agreement too, emis- 
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Line 


Uncontaminated Spectra 


Averaf 


;e Spectra 


Ref. 




HPS 


HWS 


HPS 


HWS 






(km sec -1 ) 


(km sec -1 ) 


(km sec ) 


(km sec ) 




Ha 


711 ± 10 


1606 ± 65 


676 i 19 






H/3 


762 ± 28 


1692 ± 98 


647 ± 15 


1 745 + 23 


(1) 


Br7 


702 ± 66 


1598 ± 50 


695 ± 14 


1709 ± 13 


(1) 


He I 2.06/im 


608 ± 66 


1740 ± 97 


624 ± 15 


2120 ± 74 


(1) 


PP 


619 ± 154 


1868 ± 260 


657 ± 23 


2109 ± 74 


(1) 


Pa7 






633 ± 27 




(1) 


He I 1.08/im 






646 ± 28 




(1) 


Ha 






731 ± 180 


~ 2300 ± 180 


(2) 



Table 5. Half peak separation (HPS) and half wing separation (HWS) as measured in both uncontaminated spectra and averaged spectra. 
Selection criteria for uncontaminated spectra are described in the text. IR measures are less certain due to low S/N data and h ave not 
been considered in accretion disk radii computations (see text). References: (1) this paper, (2) Gilliland, Kemper fc Suntzcff 1986j . 





R-wd 


q 


i 




ut (Xl0 10 cm) 






R ira (Xl0 10 cm) 




Fraction 




Ref. 


(M ) 


Xl0 1O cm 


«1 




Ho, 


H/3 


Br7 


1 1 CI 


H/3 


Br7 


KaicBr~f 


H/3 




0.45 ± 0.19 


0.098 ± 0.021 


0.13 


75 


1.10 ± 0.03 


0.96 ± 0.07 


1.13±0.21 


0.22 ± 0.02 


0.19 ± 0.02 


0.22 ± 0.01 


0.53 


0.48 


(1) 


1.2 ± 0.25 


0.039 


0.05 


77 


2.99 ± 0.09 


2.60 ± 0.19 


3.07±0.58 


0.59 ± 0.05 


0.53 ± 0.06 


0.59 ± 0.04 


0.93 


0.84 


(2) 


0.3 


0.055 


0.19 


75 


0.73 ± 0.02 


0.64 ± 0.05 


0.75±0.14 


0.14 ± 0.01 


0.13 ± 0.01 


0.15 ± 0.01 


0.43 


0.38 


(3) 


1.1 




0.05 


80 


2.78 ± 0.08 


2.43 ± 0.18 


2.87±0.54 


0.55 ± 0.04 


0.49 ± 0.06 


0.55 ± 0.03 


0.90 


0.81 


(4) 


0.82±0.10 




0.07 


75 


1.94±0.06 


2.04±0.08 




0.42±0.26 


0.26±0.17 




0.63 




(5) 



Table 6. The dis k parameter s determined from the line peak and line wing separations of the uncontaminated spectra. Values for q assume 
M 2 = 0.058M Q ( fcmak 1993| ). Columns 11 and 12 are the outer disk radii normalized to the primary Roche lobe radius. References: (1) 
Smak (1993), (2) Spruit & Rutten (1998), (3) Cheng et al. (1997), (4) Gilliland, Kemper and Suntzeff (1986), (5) PAPER I. Values from 
PAPER I were computed applying a different method, i.e. assuming that the outer radius extends to the 3:1 tidal resonance radius and 
deriving Mvi/Dj Rim and R otlt . The orbital inclination were taken from Smak (1993). Inner and Outer radii corresponding to the HWS 
and HPS measured on the He I and Pa/3 emission lines were not computed because of their larger uncertainties and the fact that they are 
statistically identical to the Ha and Br7 measurements. 



sion lines with larger HWS value are probably affected by 
high ve locity components from the hot-spot. HWS and HPS 



measur ;d in the averaged spectra can give information about 
the S-wave bias when compared with wing and peak sep- 
arations from uncontaminated spectra. They also test the 
assumptions used in the accretion disk radii computations. 

We compared HWS values from averaged and uncon- 
taminated spectra after correction for the white-dwarf or- 
bital motion. HWS from uncontaminated spectra, and HWS 
from averaged spectra, Kl subtracted, should be equal in 
the absence of hot-spot high velocity components. Table ^ 
shows that HWS values from uncontaminated spectra and 
Kl subtracted averaged spectra agree very well each other 
in the case of H/3 and Br7 emission lines; while Ha and He I 
2.06/^m have HWS values from Kl subtracted average spec- 
tra which are ~ 450 and 250 Km sec -1 larger than expected 
in the absence of high velocity components in the hot-spot. 
Such components bias both radial velocity curve measures 
and the accretion disk radii computation. The fact that we 
observe the largest hot-spot high velocity component in the 
Ha line is consistent with figure ^ and the less energetic 
emission forming downstream in the accretion disk. High ve- 
locity hot-spot components in the Pa/3 emission line should 
be taken cautiously, because of larger uncertainties in the 
measurement. 

It is worth noting that each average spectra was Kl sub- 
tracted by the correspondent white-dwarf orbital motion as 
determined in section [| The Kl values from radial velocity 
curve fitting are biased by the hot-spot and do not reflect 
the true orbital motion of the white-dwarf. However, be- 



cause the hot-spot b ias is believ ed t o increase the true Kl 
value (see sect ion |^, smak 197C , and Paczynski, Piotrowski 



Turski 1968 ) , then the applied correction to our averaged 



spectra should be over-estimated. 

Table ^ also reports HWS and HPS as we measured 
by hand from the average Ha spectra in figure 1 of GKS. 
Comparison of the GKS determinations with those from our 
averaged spectra gives a chance to check for changes in the 
quiescent accretion state of WZ Sge. The data in GKS were 
obtained almost one year after the last outburst of WZ Sge 
in Dec 1978 and ours were obtained about 18 years after 
the outburst. We measured larger HWS and HPS values in 
the GKS spectra than in our averaged spectra. In particu- 
lar, in the GKS spectra the HWS is ~200 Km sec -1 larger 
than in our averaged spectra; while, the HPS matches the 
values we measured on our uncontaminated spectra. There- 
fore, we may conclude that: i) the Ha hot-spot emission 
high velocity component was possibly larger one year after 
the super-outburst than 18 years later; ii) the hot-spot emis- 
sion was probably weaker at the time of GKS observations 
such that it didn't affect the accretion disk peak separation; 
Hi) the peak separation in GKS spectra may be considered 
free from hot-spot bias and implies constant outer accretion 
disk radius in the Ha emission. 

Particularly interesting is the conclusion of point Hi) 
which claims a constant accretion disk outer radius. Indeed, 
the same HPS value has been observed in all quiescent spec- 
tra of WZ Sge obtained in the last 40 years (see table ^). 
One possible explanation is in the assumption of the outer 
accretion disk radius at the 3:1 tidal resonance radius. A 
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Table 8. Results of measurements of the outer disk velocity span 
ning 40 years. Outbursts occurred in June 1946 and Dec 1978 



motion of the inner accretion disk radius toward the white- 



years 



sin ce the previous outburst. 



ri957Q~ '21 (Krzcminski fc Kraft (1964ft. C' S' Gillilnnd 
Suntzeff (1986| ), (4) fcmak (1993ft , (5) 



Kemper & 



Ncustroev (1998), (6) This 



dwarf surface (in section 4.2 we showed that the outer disk 



References: (1) Girecnstcm radius remains constant during quiescence). 



paper. 
Line 


Date 


Ayrt 


Vouter (km/sec) 


Ref 


H/3 


1956 Aug 15 


10 


720 


(1) 


H 7 


1956 Aug 15 


10 


710 


(1) 


113 


1962 Aug 4 


16 


720±14 


(2) 


Hi 


1962 Aug 4 


16 


720±14 


(2) 


Ha 


1979 July 16 


1 


731 


(3) 


H/3 


1983 Nov 24-28 


5 


730±30 


(4) 


H7 


1983 Nov 24-28 


5 


730±30 


(4) 


Ha 


1994 May 30 


16 


747±40 


(5) 


Ha: 


1996 July 27/28 


18 


711±10 


(6) 


H/3 


1996 July 27/28 


18 


762±28 


(6) 



constant accretion disk outer radius has some important im- 
plications relat ing to current outburst theory. Some d i sk out- 
burst t heor ies (Meyer-Hofmeister, Meyer & Liu 1996; Wynn 
et al 



DOOj) predict that, during quiescence, the outer ra- 
dius slowly increases up to the 3:1 resonance radius. Once 
reached at the 3:1 resonance, quiescent super-humps should 
be present. 



4.3 Smak's method of determining Rj n /R out 

The accretion disk radii ratio, K—Ki n /K ut, is often the only 
value provided by authors because of the large uncertainties 
affecting the determination of the orbital inclination and the 
binary star masses. A com monly applied method to compute 
R is the one developed by jmak (1981 ). 

Smak's method assumes axially symmetric accretion 
disk in Keplerian motion and a power law flux distribution 
(/ ~ r ~ a ) f° r the disk emission. The R and a values de- 
pend on the parameters U, As4, and A41 which are defined 
as follow: U=res/AA, with res=instrumental resolution, 
A 8 4 = log Wo. s - log Wo. 4, and A41 = log Wo. 4 — log Wo. 1, 
where Wo. 8, Wo. 4, and W0.1 are the emission line width at 
the fractions 0.8, 0.4, and 0.1 of the peak height above the 
continuum, respectively. 



Mennickent & Arenas (1998) applied Smak's method to 
WZ Sge and determined A 8 4 = 0.07, A41 =0.1, and R=0.3 
(see table ^). A ratio R=0.3 is 50% larger than the value 0.2 
we determined in section 4.1. We investigated such a differ- 



ence to understand whether it reflects changes in the quies- 
cent accretion disk size or it implies inconsistencies between 



the Smak's method and the method described in section 4.1 
We applied Smak's method to our Ha and H/3 uncontami- 
nated spectra and found Agi and A41 values matching those 
by Mennickent and Arenas (see table ^). However, we expect 
a ratio 0.2 < R < 0.3. Then, we may conclude that: i) the 
two methods are consistent and, in particular, the assump- 
tions of axisymmetric disk and power law flux distribution 
required by the Smak's method are not fundamental to the 
determination of R and Smak's method may be applied also 
to asymmetric accretion disks such as WZ Sge; ii) the two 
accretion disk radii ratio, R=0.3 by Mennickent and Are- 
nas and R=0.2 by us, may correspond to a real change in 
WZ Sge's quiescent accretion state and indicate an inward 



5 SUMMARY AND CONCLUSIONS 

We observed both hot-spot and accretion disk line emission 
to vary in shape and strength throughout an orbital period, 
and conclude that both the hot-spot and the accretion disk 
are asymmetric and anisotropically emitting. We determined 
the gas at the impact region to have just the stream velocity, 
and the accretion disk is a low density gas with little drag 
effects on the in-falling material from the secondary star. 
There are different Balmer decrements at the impact region 
and in the rest of the accretion disk, and we conclude that 
the gas has a different opacity, r, within the two regions. In 
particular, the hot-spot is optically thick in the lines, while 
the accretion disk appears to be optically thin in the lines. 

Our measured radial velocity curves, from five emission 
lines in the optical and in the infrared are found to present 
an inconsistent set of system parameters. We explain this 
result as a bias due to the hot-spot emission. In particular, 
we showed that the hot-spot emission delays the apparent 
time of the secondary inferior conjunction with respect to 
photometric phase zero. We also show that the phase off- 
set depends both on the excitation potential energy of the 
considered emission line and on the optical depth of the gas 
in the accretion disk. Our extended analysis of a selected 
sample of dwarf-novae shows increasing phase offset with 
shorter orbital period, consistent with the idea of optically 
thin accretion disks in low mass transfer rate systems. 

Summarizing our results from the emission line profile 
analysis and radial velocity curve computations presented 
here, we can formulate a probable structure for the WZ Sge 
accretion disk. The accretion disk is generally a low density, 
optically thin gas. The accretion disk is neither symmetric, 
nor has it a uniform gas density or temperature structure 
(see also PAPER I). The hot-spot is optically thick in the 
emission lines, not symmetric in shape, and not an isotropic 
emitter. The hot-spot emission arises from a multicompo- 
nent extended region in which each emission component is 
not visible at all phases and the hot-spot emission varies in 
strength throughout the orbital period. These variable hot- 
spot emissions lead to an extended hot-spot region which 
does not emit strongly when viewed from the down-stream 
direction (phase 0.5 and following), and shows a temperature 
gradient along the stream trajectory. Figure ^ summarizes 
our conclusion. 

The qualitative picture of WZ Sge provided by our ob- 
servations and data analysis shows a peculiar quiescent ac- 
cretion disk dissimilar to others observed at present. WZ Sge 
is the first short orbital period dwarf-nova showing evi- 



dence of optica lly thi n emission both in the continuum ( Cia- 
and Skidmore 1998), and the lines. Multi- 



rdi et al. 1998 



color photometry of OY Car, Z Cha, and HT C as (Wood 
~ [Wood, 199p| , |Marsh 1987] ) shows ev- 



Horne & Vennes 1992 



idence of quiescent accretion disk optically thin in the con- 
tinuum but thick in the emission lines. This work estab- 
lishes a first step toward an understanding of the spatial flux 
and material distribution within the WZ Sge accretion disk 
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Investigating the structure of the accretion disk in WZ Sge. 
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Table 9. Our measurements of the disk parameters from the uncontaminated spectra using the method described in Smak (1981) compared 
to the measure ments of Mcnnickent fe Aren as (1998). Measureme nts made assumin g a Keplerian disk f , and the data of Neustroev (1998). 
References: (1) Mcnnickent & Arenas (1995), (2) This paper, (3) |^custroev (1998} ). 



Line 


A 8i 


An 


U 


a 


1! 


Ref. 


Ha 


0.07±0.02 


0.104±0.001 


0.12 


0.5 


0.30 


(1) 


Ha 


0.072 


0.113 


0.013 


>1 


>0.2 


(2) 


H ' 


0.077 


0.122 


0.016 


>0.75 


>0.23 


(2) 


Ha 










0.20 


(2)t 


H • 










0.20 


(2)t 


Ha 








1.66±0.19 


0.08±0.04 


(3) 



and ot her TOAD (Tremendous Outburst A mplitude Dwarf- 



Howell, Rappaport & Politano 1997) c and idates. 
The accretion disk radii determined in sec. 4.1 show that 
determination of the inner and outer radius strongly depend 
on the assumed white-dwarf mass and the binary orbital in- 
clination. Unfortunately, white-dwarf mass and radius de- 
terminations are still affected by large uncertainties. The 
uncertainties on the WZ Sge white-dwarf parameters pre- 
vent us from uniquely determing the actual size of the ac- 
cretion disk, the fraction of the primary Roche lobe it fills, or 
its exact shape. Whether the accretion disk extends down 
to the white-dwarf surface or is a ring- like accretion disk 
as sug gested by some previous studies (iMever-Hofmeister 



Meyer , jz Liu 1998[ [Mennickent fe Arenas 1998; , |Howell et al 
199s(i 



lannot be uniquely determined. We determined the 
accretion disk radii ratio in WZ Sge, R=Ri„/R ollt , to be 0.2, 
and observed that the outer accretion disk radius does not 
vary during the quiescence period. One explanation for the 
constant accretion disk outer radius is the assumption that 
the accretion disk extends to the 3:1 tidal resonance radius. 
One possible test of this idea would be to search for and find 
quiescent super-humps. 
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